Thermal Equilibrium as an Initial State for Quantum Computation by NMR 
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We present a method of using a nuclear magnetic resonance computer to solve the Deutsch-Jozsa 
problem in which: (1) the number of molecules in the NMR sample is irrelevant to the number of 
qubits available to an NMR quantum computer, and (2) the initial state is chosen to be the state 
of thermal equilibrium, thereby avoiding the preparation of pseudopure states and the resulting 
exponential loss of signal as the number of qubits increases. The algorithm is described along 
with its experimental implementation using four active qubits. As expected, measured spectra 
demonstrate a clear distinction between constant and balanced functions. 



INTRODUCTION 

Liquid state nuclear magnetic resonance (NMR) [ij has 
been used to demonstrate basic concepts of quantum in- 
formation processing, including the realization of known 
quantum algorithms on a small number of qubits [2-7]. 
NMR has been leading the field both in terms of num- 
ber of qubits and in terms of coherent control techniques. 
Never-the-less, liquid state NMR quantum computing is 
widely considered to be a futile effort because the prepa- 
ration of certain initial states destroys an exponential 
amount of signal [8] and because it is not clear how to 
extract useful information by measuring expectation val- 
ues. 

Here, we present concepts that make it possible to di- 
rectly use the high-temperature thermal equilibrium den- 
sity operator as an initial state for quantum computing. 
We also show how to nnitarily evolve this thermal state 
into a state from which useful information can be ex- 
tracted. Clearly, a key ingredient of quantum informa- 
tion processing is the availability of an accessible (and 
controllable) state space that grows exponentially with 
the number of qubits and the ability to create superpo- 
sitions of these states. This requirement applies both to 
pure and mixed state quantum computing. In the case 
of pure states, the set of all superpositions necessarily in- 
cludes entangled states. However, in mixed state quan- 
tum ensembles, the ability to create (coherent) superpo- 
sitions of an exponentially increasing number of states 
does not necessarily imply that the density operator de- 
scribing the ensemble must be non-separable fs^ at each 
point in time. In fact, we argue that the concept of sep- 
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arability is of little relevance (if any) for the question 
whether a given experimental implementation of a quan- 
tum algorithm is quantum or not. Furthermore, it is 
important to note that for ensemble quantum computing 
algorithms which are based on coherent superpositions of 
basis states, the accessible exponential state space is not 
limited by the number of molecules in the sample. 

Here, we discuss a version of the Deutsch-Jozsa algo- 
rithm [10-12] for an NMR quantum computer that does 
not try to directly mimic a pure state algorithm. This 
demonstrates that a new set of algorithms is needed for 
this kind of computer. So far, there has not been an easy 
way to describe algorithms or unitary evolution on the 
density matrix of mixed states that is meaningful from 
the quantum computing point of view. Here we intro- 
duce the concept of acting on outer products, where the 
use of an extra control spin allows us to act only on one 
side of the outer product. This allows us to preserve and 
use phase information. We show that the signal-to-noise 
ratio of the presented algorithm does not diminish as the 
number of spins per molecule is increased, provided the 
total number of molecules in the ensemble is constant. 
We show that this approach provides an accessible state 
space that is of exponential size in the number of spins 
per molecule. In addition to the theoretical discussion, 
we present the first experimental implementation of this 
version of the Deutsch-Jozsa algorithm that does not suf- 
fer from the scaling problems associated with the prepa- 
ration of pseudopure states Q . 



COMPUTATION ON ONE VS MANY 
MOLECULES 

By "molecule" we mean an isolated quantum system. 
In the first subsection we address computation on one 
molecule that is restricted to be in basis states, this will 
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correspond to classical computation. Next, we address 
computation still on a single molecule that is allowed 
to be in superposition states of the chosen basis states, 
this will correspond to pure state quantum computa- 
tion. After this, we will talk about computation on many 
molecules in the two cases covered for the single molecule 
case. We will see that NMR is capable of quantum par- 
allelism over an exponential state space that is allowed 
to be in coherent superposition. 



Computation on many molecules 

In this subsection we turn to computation using many 
copies of the same molecule. The density matrix for each 
case is simply defined as the average of the density ma- 
trices of the individual molecules in the ensemble. 

Computation on diagonal density matrices 



Computation on a single molecule 

Computation on a single molecule restricted to basis states 

In this case we have a state function that starts in 
some basis state (from the exponential number of com- 
putational basis states) say state |0). A computation 
consists of logical gates which correspond to transitions 
from the initial basis state to some other basis state 
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via permutations only, where N is the size of the Hilbert 
space. The spins are never allowed to be in a superpo- 
sition state. The initial density matrix = li)OI 
for this system is a diagonal matrix with a single non- 
zero entry. During the computation, the single non-zero 
entry is permuted from one place to another on the di- 
agonal. This corresponds to classical computation and is 
not more powerful than a classical computer even though 
a quantum system is used. 



Computation on one molecule in superposition of basis states 

Here we have a system that is allowed to be in super- 
position of the form 
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Acting with a quantum gate (general unitary transfor- 
mation) on a superposition is referred to as quantum 
parallelism where, in general, an exponential number of 
computational threads take place in parallel. The den- 
sity matrix of this system is of the form where 
non-diagonal elements are non-zero in the presence of su- 
perpositions of the basis states. The control and manip- 
ulation of superpositions is a key ingredient that makes 
quantum computers more powerful than classical com- 
puters. This type of computation is referred to as pure 
state quantum computing in the literature [isj . 



Here we consider the case where the initial state of each 
molecule in the ensemble is a basis state. If the initial 
state of each molecule happens to be the same then we 
have multiple copies of the case discussed above and the 
ensemble behaves in exactly the same way. In general 
however, each molecule may be in a different initial basis 
state and the density matrix of the ensemble is diagonal 
and may have more than one non-zero diagonal entry. 
Consider an example with two molecules with one spin 
each with the following state functions 
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The corresponding density operator is 
1/2 
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If the ensemble is operated upon by permutations 
which keep the density matrix diagonal, this will cor- 
responds to parallel classical computation, where each 
molecule may be considered as a separate classical com- 
puter 
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Note that the number of molecules is very important 
in this case since this is where the parallelism is coming 
from. If each molecule contains n spins, the dimension 
of the density matrix is 2" x 2" and in order for the 
diagonal of the density matrix to be fully occupied, at 
least 2" molecules will be required [3]. 

The same conclusion holds for any ensemble that is 
described by a diagonal density operator throughout the 
computation. A diagonal density operator does not re- 
quire that all molecules are in basis states, as off-diagonal 
matrix elements can average to zero. For example, for 
two states given by: 

|^«) = i=(|o) + |i)),with|V'W)(VW' ^ ^ 
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the density operator is also given by 
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Computation on off diagonal density operators 

OfF-diagonal elements of the density operator can 
only occur if a subpopulation of the molecules exists 
whose state functions are superpositions of the computa- 
tional basis states. For example, in an ensemble of two 
molecules containing one spin each, the density operator 

1/11 
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requires that for each molecule, the state function has 
the form = -l^e'^^dO) + |1)) with arbitrary phase 

factors e'"" for fc = 1,2. 

We note that an ensemble of M identical molecules 
where each molecule contains n spins has a density op- 
erator that is the same size as the density operator of 
a single molecule and behaves in the same manner un- 
der application of unitary transforms. If it is possible 
to prepare a density operator with off diagonal elements 
and manipulate these terms, one will have a computer 
as powerful as a pure state quantum computer except 
for the measurement step which we shall address below. 
In contrast to the previous example, we note that what 
is important for the size of the available state space is 
the number of spins per molecule and not the number of 
molecules in the ensemble. 



The diagonal entries of the density matrix, i.e. when 
r = s, represent populations of the quantum basis states 
ij. Non-zero off-diagonal elements of the density matrix 
represent coherent superpositions of states. 

If all the state functions If/j*^*^^) are basis states, i.e. the 
quantum system k is not in a superposition of the basis 
states, the product Cr'^-'cs^'^' must be zero. Hence a nec- 
essary (albeit not sufficient) condition for an off-diagonal 
element p^s to be non-zero is that molecules exist whose 
individual quantum systems are in a superposition of the 
basis states \r) and |s), i.e. for each of these molecules 
the state jV'^'^-') is of the form given in equation ([2]) with 
Cr*^-* 7^ and ci'^-' ^ 0. The ensemble is said to contain 
coherence [1] between the basis states |r) and \s) if the 
sum of the terms Cr'^'c^''^'' over all molecules is non-zero, 
c.f. This implies that the available state space is 

exponential in the number of spins per molecule. One 
can compute with the diagonal elements only flij, how- 
ever even in the best case when each molecule is in a 
different basis state, the size of the state space available 
for computation is bounded by the number of molecules 
in the sample [l5| . Indeed in this case we simply have 
computation by classical parallelism. This contrasts to 
computing with off-diagonal elements [161] as we shall see 
in the following. 



ENSEMBLES OF ISOLATED QUANTUM 
SYSTEMS 

In NMR implementations, the ensemble consists of M 
identical molecules. In every molecule, the nuclear spins 
form an isolated quantum system. Here, we assume that 
each molecule contains n -I- 1 spins. (In the following 
sections, we will explain why we chose n + 1 rather than 
n spins per molecule). The resonance frequency of each 
spin is LOi — — 7;-Bo where 7; is the gyromagnetic ratio of 
spin / and Bo is the strength of the magnetic field. The 
state of the spin system corresponding to molecule k is 
given by a wave function 

j=o 

where \ j) are the standard basis states used in quantum 

(k) 

computing, Cj are the corresponding amplitudes and 
N — 2"+^. The density operator of the ensemble is given 
by 

k=l 

Each matrix element of the density operator is given by 

M 

c^}^c*S^\ (2) 
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THERMAL EQUILIBRIUM 

The thermal equilibrium of an ensemble of spin sys- 
tems is described by the Boltzmann distribution where 
the probability of the system being in state \r) is given 
by 

where Er is the energy of the r'^ eigen state of the Hamil- 
tonian of the system, k is Boltzmann's constant and T is 
temperature. 

If we assume that there are no coherences at thermal 
equilibriumjlj] the density operator of the system can be 
written as [1| 

_ exp(-g/kr) _ J_ _ H_ 
^ rr(exp(-iJ/kr)) ~ iV^ kT' 

for ||iJ|| < kT. 

In a system where the size of the couplings between 
the spins is much less than the resonance frequencies w/ 
of the individual spins the thermal density operator 
can be approximated by [IJ] 
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with ai = ^ and 

Iiz ^ ^1 (E) . . . (E) 1 (E)(J^ (g) 1(g) . . . (g) 1, 

where the PauU matrix az appears as the l*^ term in the 
product. 

UNITARY EVOLUTION, MEASUREMENT AND 
INITIAL STATE PREPARATION 

AU the topics of this section are standard in the NMR 
hterature, see e.g. [1]. Apphcation of a unitary transform 
U, whether by apphcation of rf-pulses or by time evolu- 
tion of the couphng Hamiltonian, to a density operator 
is obtained from 

p' = UpUl 

The expectation value of a Hermitian operator A is 

{A)^Tr{A p). 

For example, for the thermal density operator pth, the 
expectation value of — Ym^=i ^iz is 

{Fz)th = Tr{F, pth) 

^ n+l 

= ^ TriFz -FzY^ai hz) 
1=1 

-. n+l 

1=1 

^ n+l 

1=1 

where the facts that Fz is traceless, Tr{l) — N and if^ — 
jl were used. Note that as the number n of spins per 
molecule increases, so will the magnitude of the measured 
signal {Fz)th for ai > 0. In contrast, for the density 
operator 

Po = + (4) 
the expectation value of Fz is independent of n: 

{Fz)^TriFzPo)^^, (5) 

which is identical to {Iiz) = Tr{Iiz po). (In NMR, po 
can be created from pth using standard procedures, e.g. 
by a combination of unitary transforms and pulsed field 
gradients.) 

Application of the Hadamard transform (or of a 90° 
pulse) to pq results in the density operator 

Pi = -^(1 +"1-^1^;), (6) 

for which the expectation value of F^ (and of Ii^) is also 
ai/4. 



SCALING BEHAVIOR 

As discussed in the previous section, for pi the expec- 
tation value (Iix) is independent of the number of spins 
per molecule. Hence, the signal-to-noisc ratio for a re- 
solved resonance line of a molecule containing 3 spins is 
the same as for a molecule with 10'' spins, assuming the 
same number of molecules is in the sample, i.e. if the 
molar concentration is the same. 

As we saw earlier, there are no scaling problems with 
preparation of pi as an initial state for the computation 
starting from the thermal density operator pth- Hence, 
a scalable mixed-state based quantum algorithm can be 
constructed if the following two conditions hold: 

(1) Use of pi as an initial state for the computation. 

(2) The decision about the problem being solved is 
based upon the expectation value of Iix, where Iix is 
either parallel or orthogonal to the traceless part of the 
final density operator. 

An algorithm meeting these two conditions would over- 
come the arguments against scalability of NMR quantum 
computing [8|. (However, just as in pure state quantum 
computation, still a large number of practical or techno- 
logical impediments for realization of large-scale quan- 
tum computers would remain, such as losses due to de- 
coherence etc.) In the following section we will describe 
such an algorithm. 

THE ALGORITHM 

The starting state for the algorithm is formed by the 
density matrix pi (c.f. Eq. ([6])), the traceless part of 
which is proportional to Iix ■ As indicated in the previous 
sections, pi can be prepared from the thermal density 
operator pth without any loss of signal as a function of 
the number of spins per molecule. 

The key observation that allows us to use Iix in the 
computation is that Iix is the sum of all outer products 
which differ only in the state of the first spin: 

I ^'-^ 

^1- = 2 E |0,J>(1,J| + |1,J)(0,J| 
j=o 

where N' = N/2. Note that j runs from to iV' - 1, 
hence the size of the state space available for calcula- 
tion is N' = 2" which is exponential in the number of 
spins per molecule and is independent of the number of 
molecules in the sample. We will exploit this structure 
and use it to apply unitary transforms to the outer prod- 
ucts representing the states where spin 1 is in the state 
|1) only. Application of a unitary transform U (not con- 
trolled by spin 1) to Iix results in a sum of outer products 
where the ket and the bra do not contain information 
that we can directly use for the Deutsch-Jozsa problem. 
Instead we consider the effect of using controlled unitary 
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operations cU. For a unitary transform U over n spins, 
it is not hard to construct the unitary transform cU over 
71 + 1 spins which appHes U to spins 2 to n + 1 if spin 1 is 
in state |1) and does nothing otherwise 17|. The control 
spin is unaffected in either case. 



cu\m 
cumj) 



il(^U)\l)\j) 



For the Deutsch-Jozsa problem, given the function /, 
which is promised to be either constant or balanced, we 
are given a unitary transform which acts as follows: 

Uf\j) = {-iy^''>\3). 

It is desired to find out if / is balanced or constant by a 
single application of Uf. 

The operator cUf can be used to distinguish if / is a 
balanced or constant function using pi as the starting 
state. There are two steps to the algorithm: 

1. Apply cUf to pi to obtain the density operator p2 = 
^(1 +ai cUfh^cU]). 

2. Measure Ii^, i.e. find the value of Tr{IixP2). A 
value of ai/4 will indicate that / is constant and 
0, a value of —a/4 will indicate that the function 
is constant and 1, and a value of indicates that / 
is a balanced function as we shall see next. 

For the first step, apply cUf to pi to get 



p, = — (1 



ai cUflixcUj) 



1 

N 



(1 



N'~l 

^ Y,{cUf\0,j){l,j\cU} 



cUf\l,j){0,j\cU}}) 
1 

N'' ' 2 



N'-l 

(1 + ^ ^{|0,j)(l,j|(l2®C/]) 



(l2®C//)|l,j)(0,j|}) 

N'-l 

(-iK(^'|i,j>(o,j|}). 

The second step is the measurement of Iix which pro- 
duces 



Tr{p2hx) 



4 



if / is constant and 
- ^ if / is constant and 1 
if /is balanced. 





As in the pure state case, the measurement produces the 
same answer every time it is performed. Thus, the algo- 
rithm is deterministic. 



EXPERIMENTS 

Here, an experimental demonstration of the algorithm 
is presented. The experiments are based on the spin sys- 
tem of BOC-(i3C2-^^N-2D°-glycine)-fiuoride 18] which 
forms a five-qubit system. In the present application, 
only four of the five qubits were used explicitely in the im- 
plemented algorithm (although it was necessary to elim- 
inate interactions of the fifth qubit with the qubits that 
play an active role in our experiments). For convenience, 
here we will use the following labels for the qubits: the 
nuclear spins of the carbonyl ^^C , the aliphatic ^'^C", 
the ^^F, the ^^N and the amide atoms are labeled 1, 
2, 3, 4, and 5, respectively. 

For our demonstration experiments, we chose the fol- 
lowing test functions on n = 3 qubits: The constant test 
function was 



k{x2,x:i,X4) = 
and the balanced test function was 

fbix2,X3,X4) = X2X3 ® X4- 



(7) 



(8) 



For convenience, from now on we will drop the identity 
operator from the description of the density matrix and 
set ai = 1. Starting from thermal equilibrium, the (re- 
duced) density operator po — Iiz is prepared by a se- 
quence of INEPT- type transfers and pulsed field gradi- 
ents starting from magnetization (proportional to I^z) 
as explained in Ref. [T^j. Application of a 90° pulse to 
Po yields the desired initial state pi = Iix . 

Practical pulse sequences that implement the unitary 
transformation cUf corresponding to a given function can 
be derived from an effective Hamiltonian 



n 



log cU 
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(9) 



which effects the desired transformation in time r. 
Note that there is an infinite number of solutions to 
Eq. ([9]) from which a convenient effective Hamil- 
tonian 'H'^^ can be chosen that is compatible with 
the coupling topology of the available spin system 
[H. For the balanced test function /b defined in Eq. 
([5]) , cU is a diagonal unitary matrix in the stan- 
dard computational basis and its diagonal elements are 
{1,1,1,1, 1,1,1,1, 1,-1,1,-1, 1,-1,-1,1}. Aconve- 
nient effective Hamiltonian corresponding to cUf^ is 



nf = ^{^1 - 3/1, - h. - h. - 2/4, + 2hj., 

2I1J3Z + 2l2zl3z + ^hzhz - ^hzhzhz}- 

All product operator terms in 'H'^j^ mutually commute 
and hence can be implemented independently and in ar- 
bitrary order. The term which is proportional to the 
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FIG. 1: Schematic representation of the decomposition of 
cUf^ given in Eq. (10). The symbols are defined in the inset. 
Up to local spin operations (see main text for details), the 
propagators A-E correspond to the pulse sequence elements 
A-E of the streamhned pulse sequence shown in Fig. 2. 



identity operator can be neglected because it only con- 
tributes an irrelevant overall phase factor to the trans- 
formation cU . A schematic representation of the trans- 
formations corresponding to this decomposition of t//^ is 
shown in Fig. 1. 

The z rotations corresponding to the linear terms in 
H^-^ can be implemented by composite pulses [31 or 



by phase shifts [18|. In the spin system of B0C-(^^C2- 
^^N-^D2-glycine)-fluoride, also the terms proportional 
to Iizl2z, hzhz, and l2zl?,z can be implemented in 
a straight-forward way, because the required couplings 



Ji2 = 65.2 Hz, Ji3 = 366.0 Hz, and J23 = 67.7 Hz [18| 
are non-zero and unwanted interactions can be removed 
by decoupling. However, qubit 1 and qubit 4 are not 
directly coupled, which requires the experimental simu- 
lation of the term proportional to hzliz- The present 
implementation of this bilinear term is based on the fol- 
lowing transformation of the term l2zliz for which the 
experimental coupling constant is J24 — 13.5 Hz: 



hzhz = V{hzhz)V- 



(10) 



with 



7r 7r 
V = exp{-i-/ia;}exp{-i7r/i2/22}exp{-i-(/iy -I- hx)} 



exp{-i7r/i2/2z}exp{-i-/2y}. 

Finally, the trilinear term in the effective Hamiltonian 
which is proportional to hzhzhz was synthesized based 
on the transformation [201 



2hzhzhz = W{hzhz)W~ 



(11) 



with 



TT TT 

W = exp{-i-/i^} exp{-i7r/i^/32} exp{-i-/iy}. (12) 

The resulting pulse sequence elements can be simplified 
by standard procedures [l8|, e.g. by eliminating adjacent 
180° pulses of the same phase or by transforming some 
pulses into z rotations that can be implemented by phase 
adjustments. For example, z rotations (by angle if) can 
be implemented by a corresponding negative rotation of 
the respective rotating frame of reference. In practice, 
this results in an additional phase shift (by angle —ip) of 
all following r.f. pulses that are a ppli ed to this spin and 
of the receiver phase of this spin 



18| . The theoretically 



derived pulse sequence elements for the individual terms 
in Ti^'^ were thoroughly tested experimentally on numer- 
ous initial density operator terms. A streamlined version 
of the pulse sequence implementing cUf^ is shown in Fig. 
2. 

The durations of the pulse sequences for each term in 
the effective Hamiltonian were chosen to be integer mul- 
tiples of A = l/l^i — V2\ = 81.75 /Lis, in order to ahgn the 
rotating frames of the homonuclear ^'^C spins correspond- 
ing to qubits 1 and 2. The same pulse amplitudes and 
shapes were used for hard and selective pulses as in Ref. 
[lit ■ The bell-shaped symbols labeled "e" and "g" repre- 
sent selective e-SNOB 90° pulses [2l[ (not for the usual 
270°, but for a 90° rotation with a duration of 224 /us) 
and selective Gaussian 180° pulses (with a duration 
of 250 /xs and a truncation level of 20%), respectively 
[3]. In order to avoid nonresonant effects of these se- 
lective ^^C pulses, a standard compensation scheme was 
used [2^. For example, during an e-SNOB 90° pulse ap- 
plied to qubit 1 (offset Hz), a second e-SNOB 90° pulse 
was applied at an offset frequency of -24.462 kHz, which 
cancels the nonresonant effects at the offset frequency of 
qubit 2 at the offset -12.231 kHz. During the selective 
e-SNOB 90° pulses, qubit 3 was actively decoupled by 



applying an MLEV-4 [24 



26l | expanded sequence of 



180° pulses to spin 3 (boxes labeled "M" ). During the en- 
tire experiment deuterium decoupling was applied using 
a WALTZ-16 sequence with v^f = 0.5 kHz. 

Segments A, B, C, D, and E are related to the bilinear 
or trilinear terms of the effective Hamiltonian TH-l^ ■ Here, 
all propagators are implemented up to an overall (and ir- 
relevant) phase factor. The propagator exp{— i-| Jiz/sz} 
is implemented by segment A and additional (7r)2a; and 
{T^)bx rotations (shorthand notation for 180° rotations 
applied to qubits 2 and 5) at the end of the segment. 
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FIG. 2; Streamlined pulse sequence for the implementation 
of the unitary transformation cUf^. The durations tu — 
l/(2Jfc;) and Tki/2 — l/(4Jfcj) were chosen to be integer mul- 
tiples of A = ~ = 81.75 fj,s [13. Filled rectangles 
represent 180° pulses. The bell-shaped symbols labeled "e" 
and "g" represent selective e-SNOB 90° pulses [2l|] and se- 
lective Gaussian 180° pulses [1^, respectively. During the 
selective e-SNOB 90° pulses, qubit 3 was actively decoupled 
by applying an MLEV-4 [Sj, [2^, [2^ expanded sequence of 
180° pulses to qubit 3 (boxes labeled "M"). Dotted vertical 
lines represent local z rotations. The overall duration of the 
pulse sequence is 84 ms. 



The propagator exp{i7r 11^/22/32} is implemented by seg- 
ment B and additional {tt)2x and {tt)5x rotations at the 
beginning of the segment and {Tr)2x and (7r)iz rotations at 
the end of the segment. The propagator exp{—i-| 122/32} 
is implemented by segment C and an additional {7t)2x ro- 
tations at the beginning of the segment. The propagator 
exp{— 1-1/12/22} is implemented by segment D and an ad- 
ditional {w)4x rotations at the end of the segment. The 
propagator exp{— i7r/i2/4z} is implemented by segment 
E and an additional (TT)ix rotation at the beginning of 
the segment and a (7r)22 rotations at the end of the seg- 
ment. The additional (7r)2x and {tt)5x rotations at the 
end of segment A and at the beginning of segment B 
cancel and have been eliminated in the streamlined pulse 
sequence. Similarly, the remaining additional (tt)ix and 
{Tr)2x rotations cancel. The additional {tt)iz and (7r)22 
rotations commute with the remaining propagators and 
can be shifted to the end of the sequence, where they 
have been combined with the linear terms of Ti^^ ■ 



FIG. 3: Experimental results for (A) the constant function fo 
and (B) the balanced function ft- The spectra show the sig- 
nal of the carbonyl "C spin of BOC-(^=^C2-^^N-2D^glycine)- 
fluoride 48] corresponding to qubit 1. Our NMR experiments 
were performed using a BRUKER AVANCE 400 spectrome- 
ter with five independent r.f. channels and a QXI probe (H,C- 
F,N). The lock coil was also used for deuterium decoupling 
utilizing a lock switch. The spectra with a spectral width of 
900 Hz were acquired (128 scans) after applying the corre- 
sponding unitary transformations cUf^ and cfZ/j^ to pi = Iix- 

Application of this pulse sequence to pi ~ Iix yields 
Pi = hx — > Pb — hxliz + '^Iixhzliz + '^hxhzliz — 

'^Ilxl2zld,zliz- 

None of the product operator terms in pj, gives rise to 
detectable in-phase signal, i.e. the expectation value 
{hx) = ^T{p'^,Ilx} vanishes. 

In general, a single time measurement of {Iix) is suffi- 
cient to distinguish balanced from constant functions be- 
cause {Iix) is proportional to the integrated signal of the 
detected qubit in the frequency domain. If a full free in- 
duction decay (FID) is detected, it is possible to decouple 
all other qubits from the detected qubit 1, which makes 
the signal intensity of the detected spin (after Fourier 
transformation) independent of the number of coupled 
spins in the molecule. In our present demonstration ex- 
periments, decoupling was not used during the detection 
of the FID in order to show the details of the full niul- 
tiplet structure of the signal of the detected qubit 1. In 
the special case of the coupling topology of B0C-(^^C2- 
^^N-^D2-glycine)-fluoride, the operator p'^^ cannot evolve 
into any observable signal under the free evolution op- 
erator because qubits 1 and 4 are not directly coupled. 
Hence, even if a full FID is acquired (rather than a single 
point measurement of {Iix)), no signal is expected. The 
experimental spectrum shown in Fig. 3B matches these 
theoretical predictions. For the balanced function, the 
spectrum in Fig. 3A shows four resolved inphase signals 
(a doublet of doublets due to the couplings J12 and J13) 
with the expected intensity ratio of 1:1:1:1. 

In order to demonstrate that the lack of signal in Fig. 
3B is not simply due to artefacts or e.g. relaxation losses 
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and for an initial density operator 4/12/22/32 we expect 



' 4/12/22/32 



cU 



fb 



90°(1) 



32 



4 /1^/22/32 (16) 




>lir,VVvfV<ly(*VtiHt<'V^ with an expected 1 : — 1 : — 1 : 1 multiplet (c.f. Fig. 

4D). In all reasonable match is found between 

experimental results and theoretical predictions. 



FIG. 4: Set of control spectra (128 scans) of qubit 1 after 
application of a CNOT 4,2 gate (with a duration of 37 ms) 
to p'l, (A) and application of the pulse sequence of cUf^ (c.f. 
Fig. 2) to 2Iixl4,z (B). The additional control spectra result 
if the pulse sequence corresponding to cU is applied to po — 
2Jiz73z (C) and to po = 4/1^/22/3^ (D), followed by a 90° 
read out pulse. 



during the pulse sequence that implements cUf^ (c.f. Fig. 
2), we performed a series of control experiments. Ap- 
plication of the controUed-not operation CN0T^2 to 
yields 



, CWOT42 // 
Pb — ' Pb 



hxhz+hxhz+^lixhzliz—'^lixhzl; 



32 7 



(13) 

where the underlined terms give rise to a doublet of dou- 
blet with an expected intensity ratio of — 1 : 1 : : 0. 
The corresponding experimental spectrum is shown in 
Fig. 4A. 

If cU is applied to the inititial density operator term 
2IixIaz (rather than to pi ~ Iix), we find 



2/1 r/ 



IxJ-iz 



^Iix_ + hxhz 



/l,/32 - 2 /1^/22/32 , (14) 



where the underlined terms give rise to a doublet of dou- 
blet with an expected intensity ratio of — 1 : 1 : 1 : 1. 
The corresponding experimental spectrum is shown in 
Fig. 4B. As H^^ contains only longitudinal operators, it 
commutes with all density operator terms that contain 
only longitudinal operators. Hence, such density oper- 
ator terms are invariant under the action of C//^ . For 
example, for an initial density operator 2/12/32 we ex- 
pect 



2/12/32 
with an expected 



2hzhz 



90°(1) 



2/1^/32 



(15) 



-1 : -1 : 1 : 1 muhiplet (c.f. Fig. 4C) 



CONCLUSION 



We introduced the concept of acting on outer products 
where the use of a control spin allows us to act only on one 
side of the outer product. This provides a lucid approach 
to construct and analyze a version of the Deutsch-Jozsa 
algorithm that is scalable in an NMR setting, where the 
quantum ensemble is initially at thermal equilibrium in 
the high temperature limit. The accessible state space is 
of exponential size in the number of spins per molecule 
and the available state space is not limited by the number 
of molecules in the ensemble. Most importantly, there is 
no exponential signal loss with the number of qubits. 

The resulting algorithm has been implemented experi- 
mentally and a good match was found between theory 
and experiment. A crucial point in the implementa- 
tion of the unitary transformation cUf^ was the simu- 
lation of trilinear and bilinear coupling terms. The time- 
optimal simulation of bilinear [27] and trilinear [2^ cou- 
pling terms can be used to further reduce the duration 
of pulse sequence elements used to implement the algo- 
rithm. The size of experimentally accessible liquid-state 
NMR quantum computers is not limited by fundamental 
scaling problems related to the preparation of an initial 
state, but by practical obstacles, such as the chemical 
synthesis of suitable molecules, the control of experimen- 
tal imperfections and losses through decoherence and dis- 
sipation. However, similar questions are faced by all ap- 
proaches to experimental quantum information process- 
ing. The fact that the Deutsch-Jozsa problem can be 
decided by a probabilistic classical algorithm should not 
be taken as a point against NMR. Since NMR technology 
can be used to decide on this problem deterministically 
by a single evaluation of the function, this should be an 
encouraging point to further study NMR and its capacity 
to solve other problems. It is an interesting open ques- 
tion if in general the thermal state can be used as an 
initial state for other algorithms and whether measuring 
expectation values of spin operators is useful for solving 
other problems. 
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